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ABSTRACT 

We investigate how the shapes and angular momenta of galaxy and group mass dark matter halos 
in a ACDM iV-body simulation are correlated internally, and how they are aligned with respect to 
the location and properties of surrounding halos. We explore these relationships down to halos of 
much lower mass (10^^ hr^ M©) than previous studies. The halos are triaxial, with c/a ratios of 
0.6 ± 0.1 and a mean two-dimensional projected ellipticity of (e) = 0.24. More massive halos are 
more flattened. The axis ratios rise out to 0.6 r^-„, beyond which they drop. The principal axes, in 
particular the minor axes, are very well aligned within 0.6 rvir- High mass halos show particularly 
strong internal alignment. The angular momentum vectors are also reasonably well aligned except 
between the very outermost and very innermost regions of the halo. The angular momentum vectors 
tend to align with the minor axes, with a mean misalignment of ~ 25°, and lie perpendicular to 
the major and intermediate axes. The properties of a halo at 0.4 rvir are quite characteristic of the 
properties at most other radii within the halo. There is a very strong tendency for the minor axes of 
halos to lie perpendicular to large scale filaments, and a much weaker tendency for the major axes to 
lie along the filaments. This alignment extends to much larger separations for group and cluster mass 
halos than for galaxy mass halos. As a consequence, the intrinsic alignments of galaxies are likely 
weaker than previous predictions, which were based on the shapes of cluster mass halos. The angular 
momenta of the highest concentration halos tend to point toward other halos. The angular momenta 
of galaxy mass halos point parallel to filaments, while those of group and cluster mass halos show a 
very strong tendency to point perpendicular to the filaments. This suggests that group and cluster 
mass halos acquire most of their angular momentum from major mergers along filaments, while the 
accretion history of mass and angular momentum onto galaxy mass halos has been smoother. 

Subject headings: galaxies: formation — galaxies: halos — galaxies: structure — dark matter — 
galaxies: clusters: general — methods: N-body simulations 



1. INTRODUCTION 

The three-dimensional structure of the dark matter 
halos that host galaxies, groups, and clusters is an im- 
portant aspect of their nature that can provide insight 
into their formation and affect the luminous structures 
within. The orientation of the halo shapes and angular 
momenta, both internally and with respect to surround- 
ing halos, provide important constraints on other studies 
of galaxy formation and evolution. 

Halos formed in cosmological simulations are gener- 
ally not spherical, but have an ellipsoidal shape. There 
have been several studies of the shapes of halos in low 
resolution A'^-body simulations based on the standard 
cold dark matter (CDM) paradigm (Frcnk ct al. 1988; 
Eubinski & Carlbcrg 1991; Warren ct al. 1992 (hereafter 
rW92) : .Cole fc Lacev 1996). These studies have found 
that halos are usually triaxial, with a preference for pro- 
late figures at small radii and more oblate figures at large 
radii, and have min or-to-major axis ratios rang ing from 
0.3 to almost unity. IDubinski fc Carlberd l|199l() . in sim- 
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ulations of the formation of isolated galaxies where the 
effects of the external tidal field were prescriptively su- 
perimposed, found that the projected two-dimensional 
ellipticities peak around e — 0.5, where 

e=l-<7, (1) 

for a projected axis ratio of q. IW92I found that while this 
holds in the inner regions, the location of the peak falls to 
e — 0.25 beyond 50 kpc. Th e behaviour of the axis ratios 
with radius is controversial: IDubinski fc Carlber'^ ljl991f) 
and iW92 find that the axis ratio s increase (become 
mor e spherical) w i th rad ius, while iFrenk et alJ 1)19881) 
and ICole fc LacevI l)1996j) find that they decrease with 
radiiis. IWgl " have also studied the internal alignment 
of the ellipsoid principal axes. They found that both the 
major and minor axes of halos are extremely well aligned 
out to 40 kpc. 

More recently, several authors have performed large 
high resolution simul ations u s ing th e currently-favoured 
ACDM co smolog y fBullock "200^; 'Jine fc Suto' '2Ciol 
(herea f ter lJS02fl: Kasun fc Evrard 2004 (hereafter 
lpX)|) : lHopkins et al.ll2005l fhereafter HBB05)). iBulloci] 
(2003) finds that the c/a axis ratios arc a strong function 
of halo mass, and range from 0.55 at 10^"* /i^^ M© to 0.7 
around 10^^ M©, with a distribution that is peaked 
but has a large tail to small axis ratios. He also finds that 
the inner 30 kpc of halos are more spherical than the 
outer regions, i.e . the a xis ratios decrease with radius. On 
the other hand, IJS02I find c/a axis ratios that increase 



2 



Bailin and Steinmetz 



with radius, decrease slightly with mass, and are well fit 
by a Gaussian centred at c/a = 0.55 with a width of 0.11. 
They find that the major axes of halos are relatively well 
ahgned — typically cos^n ~ 0.8, where On is the angle 
between the major axis at small or large radius compared 
to that at an intermediate radius. The align ment of the 
middle axes is somewhat poorer, but rrsol argue that 
this is due to the inclusion of nearly prolate halos whose 
axes a re degenerate and therefore not well determined. 
IKE04I find more spherical shapes in their very large sam- 
ple of high mass halos (M > 3 x 10^^ M©), possibly 
due to the s pherical outer boundary they impose, with 
c/a « 0.65. IhBRO,'^ measure the intermediate-to-major 
axis ratio h/a of cluster mass halos and find a mean of 
0.67, w ith a mass dependence similar to that found by 
IBullock ( 20021. and axis ratios that decrease with radius. 

The shapes of dark matter halos can have impor- 
tan t observ ational consequences (for a good review, 
see ISacket3 [1999) . On galactic scales, they can af- 
fect the coherence of tidal streams. Some authors have 
claimed that the thinness of the tidal stream asso- 
ciated with the Sagittarius dwarf spheroidal indicates 
that the halo of of the Milky Way i s nearly spheri- 
cal, with c/a > 0.8 lllbata et all 120011: iJohnston et all 
120051 iMartfnez-Delgado et all l2004j) . However, more 
recent studies suggest that the material that makes 
up the stream was stripped from the satellite too re- 
cently to have had time to undergo differential pre- 
cession, wh ich severely we a kens th e constr aints on the 
halo shape llHe]mill2004albll . IHelmil lj2004br even claims 
that the stream is best fit assuming a prolate halo 
elongated perpendicular to the d isk with c/a « 0.6. 
iLaw. Johnston, fc Maiewskil ()2005D also find that the ve- 
locities of stars in the leading stream can only be fit 
with a prolate halo, but that the precession of the lead- 
ing stream with respect to the trailing stream can only 
be fit with an oblate halo. These contradictory results 
suggest that evolution of the satellite orbit or other ef- 
fects of a live Milky Way potential (as opposed to the 
static potential that has been used in all of these stud- 
ies) are important for determining the shape of the Milky 
Way halo using the Sgr stream. Other measures of the 
Milky Way ellipticity using the flaring of the gas disk at 
large radius or the anisotropy of stellar velocities suggest 
that the halo is oblate with a flattening oi c/a ~ 0.8 
((oiling & Mcrrificld 2003). The shapes of the halos of 
external galax ies can be measur ed using the flaring of 
the gas layer (lOUing fc Merrifield .2000.') . the projected 
shape of X-rav gas (jBuote et al.ll2002tl. or the kinemat- 
ics of polar ring galaxies (jSackett et al.lll994fl . These 
methods suggest that galaxy halos have a wide range of 
flattenings from c/a ~ 0.3 - 0.8. Assuming that the 
shape of the stellar halo traces that of the dark mat- 
ter halo, the stacked images of edge-on disk galaxies 
in the Sloan Digital Sky Survey are consistent with a 
mean c/a = 0.6 (Zibctti et al. 2004). A new method for 
measuring the shapes of external galaxy halos is weak 
gravitational lensing. By measuring the azimuthal vari- 
ation of the shear with respect to the position angl e 
of the visible lens gal axy (jNataraian fc Refregiedl2000j) . 
iHoekstra et al.l 1)20041) detected an average projected halo 
ellipticity of (e) = 0.33 for halos with an average mass of 
8 X 10^^ M0. This detection also implies that the ori- 
entation of the visible and dark mass in galaxies must be 



similar. On group and cluster scal es, X-ray observations 
and the Sunyaev-Zeldovich effect ()Siinvaev fc ZeldovichI 
1198(1 can be dir ectly used as a probe of halo ellipticities 
(ILee fc Sutcl l2003. 2004). 

The orientation of the angular momentum in halos has 
also been studied in cosmolo gical numerical simulations . 
Early low-resolution studies (iBarnes fc Efstathioul Il987t 
iFrenk et alJ Il988t iCole fc Lacevl I1996D gave conflicting 
results due to the difflculty of measuring the direction of 
the angular momentum with few particles. Other CDM 
studies (Dubinski 1992: W92) have found that the direc- 
tion of the angular momentum at different radii is usually 
the same, but that the distribution of alignments has a 
tail that stretches all the way to ant i- alignment. This 
result is verified by recent high-resolution ACDM sim- 
ulations (Bull ock et al.. 2001al. CDM studies have also 
found that the angular momentum is most often aligned 
with the minor axis and perpendicular to the major axis, 
although there is some scatter (Dubinski 1992; W92). 
This result has not yet been thoroughly tested in high- 
resolution ACDM simulations. 

Internal misalignment of the angular momentum 
can have a number of observational conseq uences. It 
may cause galactic warps (Ostrikcr & Bi nnevl 119891: 
[Debattista fc Sellw oodI Il999t iLopez-Corredoira et alJ 
1200 2: Bailin 2004), or manifest itself in anisotropic distri- 
butions of the orbits o f satellite galax ies ( Holmbcrg^l^J^ 
i Zaritskv et aD Il997t lAubert et all [2004: .Knebe et alJ 

Going beyond individual halos, the shapes and angular 
momenta of nearby halos can be correlated due to initial 
conditions or dynamical evolution. This subject has at- 
tracted increased interest recently due to the emergence 
of weak gravitational lensing as a method to measure 
the projected mass density in front of background galax- 
ies. Intrinsic correlations between the projected shapes 
of luminous galaxies act as spurious background signals 
in weak lensing, so predicting their magnitude is impor- 
tant. At a more fundamental level, the degree of corre- 
lation between structures can be tested against models, 
and can inform our understanding of the origin of halo 
shapes and angular momenta. 

Measurements of halo alignments and correlations 
come from two sources: cluster orientations and large 
galaxy surveys. The study of the alignment of cluster 
orientations was pioneered by Binggcli (1982), who used 
the locations of the constituent galaxies to determine 
that the major axes of clusters separated by less than 
15 Mpc tend to point toward each other. While 
some authors have not found any such alignment (e.g. 
Struble fc Peebles 1985), larger samples of both galaxies 
and clusters, along with improve d error estimates, have 
confirmed th is result (e.g. Flin Il987t iRhee fc KatgertI 
ri987tlPlioni3 rr994V Similar results are seen when us- 
ing the major axis of the brightest cluster galax y as a 
meas urement of the cluster orientation (Lambas et alJ 
'1990). Enhancements in galaxy counts along the major 
axis of brightest clus ter galaxies out to 15 h^^ Mpc have 
also been detected (lArgyres et al.l llQSd: iLambas et alJ 
119881: iMuriel fc Lambad \198^ . The cl uster potential 
is better probed by X -ray emitting gas ([SarazinI Il986t 
ILee fc Suto"2003'. '20041 While early studies using X-ray 
contours found no alignment of c lusters with the large 
scale structure (|Ulmer et al.ll989j) , more and better data 
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have confirmed that the orientation of both the substruc- 
ture and the main clus ter potential tends to point toward 
neighbouring clusters l)West et al.lll995t iChambers et al.l 
11)02). 

While the principal axes of clusters can be deter- 
mined from optical or X-ray photometry, the angular 
momentum direction is very difficult to determine. In 
disk galaxies, on the other hand, the angular momen- 
tum direction of the baryons, presumed to be perpen- 
dicular to the orientation of the disk, can be measured 
much more easily than the shape of the dark matter 
halo. Although there may be some misalignment be- 
tween t he angular momentum o f the baryons and dark 
matt er ijvan den Bosch et al.ll200 2: Sharma & StcinmcQ 
12004^ ■ spiral galaxies still provide the best targets for 
detecting angular momentum alignments and correla- 
tions. Studies with small samples of spiral galaxies (less 
than a few hundred) have generally found no correla- 
tion between the orientation o f the angular momentum 
and th e large scale structure (Cerne fc Pete rsonlll99Ct 



over the same range of separations (iBarnes fc Efstathior 
19871 iWest et al.lll99lHvan Haarlem fc van de Wevgaeri 



Han et al. 1995; Cabancla & Dickev 1999; see however 



Navarro. Abadi. fc Stcinmctz 200j|). With a sample of 



618 lenticular and disk galaxie s in the local supercluster, 
iKashikawa fc Okamural l|1992() found that while the full 
sample was consistent with an isotropic distribution of 
angular momenta, those galaxies within 2 h^^ Mpc of 
the supergalactic plane tend to have spin vectors point- 
ing in the plane, while those above or below the plane 
tend to ha ve spin vectors that po int toward or away from 
the plane. iNavarro et alJ l)2004f l also find a clear excess 
of galaxies whose angular momenta lie in the supergalac- 
tic plane. Larger samples of galaxies provide further ev- 
idence of alignments between spin and the large scale 
structure: iPen. Lee, fc Selia^ l|20fl(1l) have found that the 
spin directions of the 12,122 spiral galaxies in the TuUy 
catalog are positive ly correlated a t sepa rations less than 
3 Mpc, while iBrown et alJ l)2002f ) have measured 
intrinsic correlations between galaxy orientations at a 
range of angular separations in the 2 x 10^ galaxies of 
the SuperCOSMOS survey. 



Lin ear tidal torque theory (jPoroshkevi' , 

IT981 can be used to predict the directions of the an- 
gular momentum vecto rs and their correlations with the 
surrounding structure llPen et alJl200nHLee fc Penll2nn(l 
■ II 120011 IF 



20M iCrittenden et al.l l200li Forciani et all l2002aibfl ^ 

These studies have found that the angular momenta of 
halos tend to lie perpendicular to the large scale struc- 
ture, and that the correlation of the halo spin vectors 
with each other exists but is very weak. While tidal 
torque theory performs reasonable well at predicting 
the e volution of the magni t ude of the angular m omen- 
tum ijSugerman et al.ll2000(l . |Porciani et a l. (2002a') have 
tested the predictions of spin directions against iV-body 
simulations, and found that the spin axes of A^-body ha- 
los show significant misalignment compared to the tidal 
torque predictions, with a mean misalignment of ^ 50° 
at z = 0. Therefore, A^-body simulations that take the 
full non-linear dynamics into account are necessary. 

Some early numerical work at predicting the intrinsic 
alignments using simulations with power law or CDM 
power spectra found that the major axes of cluster mass 
halos tend to point toward other nearby clusters over 
scales of ^ 15 h^^ Mpc, and that there is a very weak ten- 
dency for the major axes to be correlated with each other 



19931 iSplinter et all Il997|) . More recent high resolu- 
tion TV-body simulations in a ACDM cosmology have 
been studied to search for alignments of the major axes 
(lOnuora fc Tho"^ l20(¥i lkE04.: HRB0,5>1. the angu- 
lar momentum vect ors l | Hatton fc Nininll2001j) . or both 
()Faltenbacher et alJ 120021 hereafter |F02|) . These stud- 
ies have found a strong tendency for the major axes of 
cluster mass halos to point toward other clusters out to 
several tens (or even hundreds) of Mpc and to corre- 
late with each other out to 20 Mp c. The situation 
for th e angular momentum is murkier. iHatton fc NiniiJ 
()2001l) have found that the angular momenta of halos 
tend to lie parallel to the large scale structure, while 'FO^ 
have found that they lie perpendicular to the large scale 
structure, i n agreement with the prediction from the lin- 
ear theory ijLee fc Penll2001[) . Correlations of halo angu- 
lar momentum vectors with each other appear weak at 
best. Two groups have searched specifically for correla- 
tions between the ellipticities of z « 1 galaxies in the 
Virgo Consortium ACDM model f Jen kins et all 1199^ 
in order to predict th e effect on weak lensing signals. 
iCroft fc Metzled 1)2000(1 assumed that the shapes of the 
galaxies were identical to the shapes of the dark mat- 
ter halos, while Heavens ct al. (2000) (hereafter HRHOOf) 
made separate predictions about ellipticals, which were 
assumed to share the shape of their dark matter halos, 
and spirals, which were assumed to lie orthogonal to the 
halo angular momentum. In both cases, the correlations 
are small, but detectable. 

In this paper, we present an extensive study of the 
shapes of galaxy and group-mass halos in a large high 
resolution ACDM A^-body simulation. We study the in- 
ternal alignments of all of the principal axes and the an- 
gular momentum. We also study the alignment of all of 
these quantities with the large scale structure, and how 
they are correlated in halos of different separations. This 
work improves upon earlier studies of the internal struc- 
ture of halos by using large high resolution simulations 
in a currently-favoured ACDM cosmology, by studying 
both the alignment of the angular momentum and the 
shape, and by using a method that allows us to quantify 
our errors and therefore feel confident about the source 
of any measured misalignments. Previous external align- 
ment studies have all been restricted to massive clusters; 
we improve upon this significantly by reaching down to 
galaxy mass halos, by studying the mass dependence of 
the correlations, by studying the alignments of both the 
principal axes and the angular momenta, and by includ- 
ing the oft-neglected intermediate and minor axes. The 
structure of the paper is as follows. In §|5|we present the 
details of the simulation and describe the method used 
to measure the principal axes and angular momentum 
vectors of the halos. § 01 presents the overall shapes of 
the halos and how they change with radius. We discuss 
the internal alignment of the principal axes and angular 
momenta in § 0] while we explore the alignment of these 
quantities with external halos in § [S] We discuss what 
these results imply for the origin of halo shapes, angular 
momentum, and warped galaxies in § and summarize 
the results in §[3 

2. METHODOLOGY 
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2.1. The simulation 

The simulation used here is the same as the one used 
in'Baihn fc Steinmetd (|2004). It consists of 512^ iV-body 
particles in a periodic box of length 50 h^^ Mpc, in a 
low density flat universe (fi = 0.3, J^a — 0.7, h — 0.7, 
as = 0.9). The particle mass is 7.757 x 10^ M©, 
and the force softening length is 5 h^^ kpc. Halos were 
found using the standard friends-of-friends (FOF) algo- 
rithm with a linking length of 0.2 times the mean inter- 
particle separation. 

In order to accurately measure the direction of the 
principal sufficiently large number of particles per 

halo are required. The direction of an axis can generally 
be determined to within an angle of 

1 /r 

9cTi = — 1= radians, (2) 

where N is the number of particles used and r is the 
relevant axis ratio: b/a for the major axis, c/b for the 
min or axis, and max(6/a,c/b ) for the intermediate axis 
fsee lBailin fc St einmetz"2004') . For the purposes of mea- 
suring internal alignments, we would like the angular er- 
rors to be less than 10°. This requires on the order of 
200 particles. Since each halo is split up into 6 radial 
shells, and this accuracy is required in each shell, each 
halo should have over 1200 particles. For convenience, 
we adopt a cutoff of 10^^ Mq for the mass of the 
halo, or 1289 particles. 

There are 3869 halos in the sample with masses extend- 
ing from 10" /i-i M0 to 2.8 x lO^^ h-^ Mq. 451 of the 
halos have masses in the range 10^^ - 10^'^ M©, while 
62 of the halos have masses greater than 10^^ M©. 

2.2. Measuring the axes 

A standard technique to measure halo triaxiality in 
simulations is to use an iterative approach, where the 
particles are initially chosen to lie in a sphere or spheri- 
cal shell, an ellipsoid is fit to these particles, and particles 
are chosen for the next iteration based on the new ellip- 
soid (e.g. W92:). While this works for simulations that 
have sufficiently low resolution that overmerging erases 
substructure, we find in agreement with lJSO'^ that in high 
resolution simulations, the presence of substructure pre- 
vents this technique from converging for a large fraction 
of halos. 

JS02 have adopted a novel approach aimed at directly 
measuring isodensity contours. They assign SPH-like 
densities to halo particles, and then measure the prin- 
cipal axes and axis ratios of particles with densities near 
the nominal density of the isodensity contour. Due to the 
presence of substructure, this procedure often picks out 
disconnected shells in addition to the particles that define 
the isodensity contour of the smooth distribution. IJS02I 
use the FOF algorithm to select the largest structure 
that fulfils the density criterion and therefore eliminate 
the substructure. 

While this algorithm works well for very high- 
resolution halos, such as those I.TS02I use to demon- 
strate the technique (all of which contain iV > 6 x 10^ 
(ping & Suto 2000)), we encountered difficulties using it 
on more moderate resolution halos. In particular, wc find 
that it is not possible to find an optimal FOF linking 
length for eliminating substructure; if the linking length 



is too large, many "contours" contain obvious discon- 
nected substructures, while reducing the linking length 
sufficiently to eliminate this problem reduces it to the 
point where in many halos, the single ellipsoid corre- 
sponding to the smooth distribution is brok en up by the 
algorithm into several disconnected pieces. lJS02f s algo- 
rithm also uses a relatively small number of particles to 
determine the shape. The error in the determination of 
the direction of the principal axes of a particle distribu- 
tion goes as 7V~^/^ (see eq. |2]). Measuring the internal 
alignment of a halo requires small well-understood errors, 
and therefore as many particles as possible. 

We take the f ollowi ng approach, similar to that taken 
bv lFrenk et al.l l|1988f) . The center of mass is calculated 
iteratively in spheres centred on the center of mass of the 
sphere in the previous iteration, starting with a sphere 
containing all of the particles. The radius of each suc- 
cessive sphere is reduced by 90%, and the procedure is 
iterated 25 times, by which point it has converged. The 
particles of each halo are transformed into this center of 
mass frame, and the velocities are transformed into the 
center of velocity frame. Each halo is then divided up 
into six concentric spherical shells with outer radii R of 
1.0, 0.6, 0.4, 0.25, 0.12, and 0.06 times the virial radius 
Tvir. These radii are chosen to allow easy comparison 
to the isodensity contours of |JS02 > The outer radius of 
each shell also forms the inner radius of the next larger 
shell, except for the innermost "shell" , which is actu- 
ally a sphere extending to the halo center. One would 
like to use the inertia tensor to measure the principal 
axes of the mass distribution. However, the inertia ten- 
sor can be dominated by substructure in the outer part 
of the shell. Therefore, we weight particles by 1/r^ so 
that each mass unit contributes equally regardless of ra- 
dius (|Gerha rd'1983'). Within each shcU, we calculate this 
reduced moment of inertia tensor 

which we then diagonalize. The axis ratios a, 5, and c 
are the square roots of the eigenvalues {a > b > c), and 
the eigenvectors give the directions of the principal axes. 
There are no particles in common between radial shells, 
so the measurements of the axes at different radii are 
completely independent. We also calculate the angular 
momentum of the particles in each shell. 

To calculate the error in the axis and angular momen- 
tum determinations, we perform a bootstrap analys is of 
the particles in each radial shell ijHevl et alJll994j) . If 
the shell contains N particles, we resample the shell by 
randomly selecting N particles from that set allowing for 
duplication and determine the axes and angular momen- 
tum from this bootstrap set. We do this 100 times for 
each radial shell. The dispersions of these estimates of 
the axis ratios, directions of the axes, magnitude of the 
angular momentum, and direction of the angular mo- 
mentum around the measured values are taken formally 
as the "Ict" error of each of these quantities. 

Using spherical shells, rather than the ellipsoids de- 
fined by the isodensity contours as in IJSOI does 
not affect the orientation of the principal axes. It 
does, however, bias the derived axis ratios toward 
spherical. To calibrate the magnitude of this bias, 
we have constructed prolate Poisson-sampled NFW 
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ratio, there are 12 independent points, representing the results for 0.4- » lo'^ - 10" h-' Me 

each of six radial bins and two different density profiles, but they : + - 

are virtually indistinguishable. The solid line is the empirical fit of 0.1 

equation jlj. '/'•>• 




0.0 0.2 0.4 0.6 0.8 i.O 

c/a 

Fig. 2. — Histogram of the minor-to-major c/a axis ratio for each 
halo in the simulation, as measured in the R = 0.4 rvir shell. 



(jNavarro. Frenk. fc Whit3ll996j) and singular isothermal 
halos with 10^ particles each that have known c/a axis 
ratios ranging from 0.3 - 1. The NFW halos had con- 
centration parameters Cvir = 10. In Figure ^ we plot 
the input c/a axis ratio for each combination of radial 
bin and density distribution as a function of the c/a ra- 
tio measured for the halos constructed using the method 
described above. The points for different radii and den- 
sity distributions are virtually indistinguishable; the only 
important parameter is the input axis ratio. The solid 
line is an empirical fit to these points: 

(c/a)t,uc = (c/a),'f ,,„red- (4) 

The relationship for b/a is identical. The corrected axis 
ratios (&/a)true and (c/a)truo are used throughout the 
remainder of this paper. 

3. SHAPES 

Dark matter halos are well approximated by ellipsoids 
(jJS02) and are well described by the intermediate-to- 
major and minor-to-major axis ratios b/a and c/a. Fig- 
ure|21is a histogram of the c/a axis ratios as measured in 
the R = 0.4 Tvir shell for all of the halos in our sample. 
Unlike |jS02l we find that they are not quite Gaussian, 
but rather have a tail toward very flattened halos as seen 



Fig. 3. — Median of the b/a (top), c/a (middle), and c/b (bot- 
tom) axis ratios for each radial shell. Crosses, asterisks, and di- 
amonds represent the mass ranges lO'^^ — lO'^^ Mq, 10^'^ — 
10^^ h-^ Mq, and lO^^ - 3 X lO" h'^ Mq respectively. Error 
bars represent the la width of the distribution. The error in the 
median is typically 0.002 for the crosses, 0.005 for the asterisks, 
and 0.01 for the diamonds — smaller than the sizes of the sym- 
bols. Crosses and diamonds are offset in radius for clarity. 



by'KE04', although the tail is not as extreme as that seen 
by Bullock ( 2002) . The distribution of b/a and c/a values 
measured at other radii have a similar shape. The lack 
of very flattened halos in IJS02I may be a result of their 
exclusion of halos deemed to be interacting. The axis ra- 
tios we find are intermediate between the quite flattened 
halos found by JS02 and the more spherical halos found 
by KEM 

Early studies suggested that the coldness of the Sgr 
stream indicates that t he dark matter halo of the Milky 
Way has c/a > 0.8 (Ibata ct al ] l 2001l iJohnston erall 
2005; Martfnez-Dclgado et al. 20(3). Only 5% of the ha- 
los shown in Figure|2|have axis ratios so large. While this 
is not negligible, it is uncomfortably small, and forces us 
into the anti-Copernican situation of living in an excep- 
tional galaxy. However, new models of the Sgr stream 
that are more careful about matching observations of 
the body of the dwarf find that the stars that constitute 
the stream were stripped from the body of the satel- 
lite too recently to have had time to undergo differential 
procession, thereby severely weakening the constraints 
on the halo ellipticity ( Helmi 2004a b); current models 
are unable to simultaneously fit the velocities of stars in 
the leading stream and t he orbital planes of the lead- 
ing and trailing streams ijLaw et alJ l2005|) . It should 
also be noted that our simulations do not take into ac- 
count the effects of baryonic physics. There is some evi- 
dence tha t baryon cooling leads to more spherical halos 
(|Dubinskilll99l iKazantzidis et al.ll2004D . Therefore sys- 
tems in which most of the baryons have cooled, such as 
disk galaxies like the Milky Way, may have dark mat- 
ter shapes that are more spherical than those presented 
here. Observations of external galaxies using a variety of 
methods find halo fiattenings t hat range from 0.3 to 0.8, 
in agreement with our results l)Sackettill999(l . 
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Fig. 4. — Histogram of the triaxiality T of all halos. The solid, 
dashed, and dotted histograms represent the halos measured at 
0.12, 0.6, and 1.0 rvir respectively. 



The radial dependence of the axis ratios is shown in 
Figure El There are three distinct regions of the halo. 
Over most of the halo, the axis ratios increase with radius 
(i.e. the halos become more spherical). Near the virial ra- 
dius, infalling unvirialized structure causes the axis ratios 
to drop. In the central 6% of the virial radius, the axis 
ratios rise. However , this is probably an artif act of the 
numerical softening l|Bailin fc Steinmetd )2004'l. Further 
evidence that this is a numerical effect comes from exam- 
ining the location of this increase in sphericity for halos 
of different mass. The increase occurs at a larger fraction 
of the virial radius for low mass halos, i.e. at a similar 
physical radius. This complicated and non-monotonic 
radial dependence may explain the discrepancy between 
studies that have found t hat flattening increases with 
radiu s l|FreTik et a.lJ 11083 ICole T.a.cevI ITflM 'Bullqfia 
l2002j) and those that have found t hat it decre ases with 
radius mubin ski fc CarlberdllMt IW9l IJSnl IHBBOl . 
The mass dependence of the axis ratios is shown by the 
different symbols in Figure 01 The highest mass halos 
have smaller axis ratios at all radii than the smaller ha- 
los (jW92; Bullock 2002: JS02: KE04). Our halos extend 
down to masses an order of magnitude smaller than any 
other study; the difference between the most massive ha- 
los (the diamonds) and the galaxy mass halos is more 
pronounced than the difference between different masses 
of galaxy halos. The high mass halos show particularly 
strong flattening near the virial radius. 

As seen in Figure (31 the c/h ratio falls steadily 
with radius, indicating a transition from prolate fig- 
ures in the center to ob l ate fi gures at large radii 
ijDubinski fc Carlberd IT99ll IW9^ . Triaxiality can be 
quantified by the parameter T: 



0.0 I 

0.0 



c/a 



T 



(5) 



ijFranx et al.lll99ll) . Purely prolate halos have T — 1 
while purely oblate halos have T — Q. We have mea- 
sured T at three radii: i? = 1.0 Tvir, where infalling ma- 
terial results in substantial flattening, R = 0.6 rvir, where 
the halos are at their least flattened, and i? = 0.12 rvir, 
where the interior of the halos are at their most flattened. 
Histograms of T at these radii are shown in Figure 01 
The interior regions of halos clearly tend to be prolate 
(solid histogram). As the flattening decreases at larger 



Fig. 5. — Intermediate axis ratio 6/a as a function of minor axis 
ratio c/a for all of the halos, measured at R = 0.4 rvir- The inner 
and outer contours enclose 68% and 90% of the halos respectively. 
The thick lines have constant values of the triaxiality parameter 
T. The separations between the prolate, triaxial, and oblate pop- 
ulations occur at T = 1/3 and 2/3. 



radii, many of the halos become more oblate (dashed 
histogram), although still more are prolate than oblate. 
Near the virial radius, there is a small shift back toward 
prolate shapes (dotted histogram) . Figure O shows the 
full relationship between b/a and c/a for all of the ha- 
los in our sample, measured at R = 0.4 rvir, where the 
values are most typical for the halo as a whole. The pre- 
ponderance of prolate and triaxial halos over oblate halos 
is clearly seen. 

IHBB05I examined the evolution of cUipticity of cluster 
mass halos at a wide range of redshifts in a ACDM simu- 
lation, and found that halos of a given mass have higher 
ellipticity at higher redshift. In particular, they find that 
(e) = 0.33 + 0.05z for all of their halos, with slightly 
stronger evolution (d (e) /dz — 0.06-0.07) for their lowest 
mass halos, which overlap with our high mass halos (note 
that these authors define e = 1 — h/a). There are two pos- 
sible interpretations for this result: either individual ha- 
los become more spherical with time, or halos that form 
later are intrinsically more spherical and the evolution 
seen by HBB05 is due to the growth of individual halos, 
resulting in a typically later formation time of halos of a 
given mass as time pas ses. It is interestin g to compare 
this with the results of iBailin fc Steinmet z (2004), who 
studied the evolution of the ellipticity of individual halos 
at z « 0. Bailin fc Steinmetz (2004) find that individual 
halos have a mean d{h/a)/dt = 0.007 Gyr~^, cor- 
responding to d (e) /dz = 0.07, in good agreement with 
HBBO.^ We conclude that the shapes of individual halos 
become more spherical with time. 

Weak lensing measurements have recently begun to 
probe the two-dimensional projected ellipticity of the 
lensing mass distribution l|Hoekstra et al.l 120041) . For 
each halo in our sample, we have calculated the pro- 
jected axis ratio q using the method of iStarkI fl977*) for 
two random orientations, and calculated the ellipticity e 
using equation A histogram of the results for the 
i? = 0.4 rvir shell is shown in Figure [HI The mean and 
median of the distribution of ellipticities are Cmcan = 0.24 
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Fig. 6. — Projected ellipticity e of each halo seen from two 
random orientations, measured at i? = 0.4 rvir- 



TABLE 1 

Number of halos with axes determined to within 

0.2 RADIANS 



Radius 


Major axis 


Intermediate axis 


Minor axis 








0.06 


904 


388 


621 


0.12 


2370 


1005 


1628 


0.25 


2942 


1762 


2540 


0.4 


2545 


1488 


2421 


0.6 


2322 


1363 


2324 


1.0 


2973 


2231 


2944 



and eincdian = 0.23 respectively, with a la width of 0.13. 
This is consistent with the lower limit of (e) = 0.33^oog 
found bv lHoekstra et all l)2004D from stacked weak lens- 
ing measurements around galaxies. It is smaller than the 
elHpticities found by Croft & Metzler ( 2000) (note that 
they quantify ellipticity as e = (1 — (7^)/(l + g^), in which 
units our mean ellipticity is e = 0.27). However, most of 
the halos used in their study have higher masses than the 
galactic-mass halos studied here, and they included rela- 
tively poorly-resolved halos which may skew the results 
toward higher ellipticities. 

4. INTERNAL ALIGNMENT 
4.1. Principal axes 

We compare the alignment of the principal axes within 
each halo to see whether the approximati on of t he halo as 
a set of concentric ellipsoids is justified l|JS02D . In order 
to determine whether the axes are aligned, the directions 
of the axes must be well determined, otherwise apparent 
misalignments may be produced owing to measurement 
errors. Therefore, we restrict ourselves in this section to 
axes whose bootstrap error is less than 0.2 radians. The 
number of halos satisfying this criterion at each radius 
for each axis is given in Table 

To understand the effect of the error on the determina- 
tion of the alignment, consider two axes which are intrin- 
sically perfectly aligned, but are each measured with an 
error of 0.2 radians (note that this is the worst possible 
case — the median error of the sample is 0.1 radians). 
Due to the measurement error, these axes will appear 
to be misaligned by an angle ^spurious- The component 
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Fig. 7. — Histograms of the direction cosine between the major 
axis of the halo at i? = 0.4 rvir E'nd the major axis at = 1.0, 0.6, 
0.25, 0.12, and 0.06 rvir (top to bottom). Due to the symmetry 
of the axes, this is always positive. Each histogram contains all 
halos where the major axes at both radii are determined to within 
0.2 radians. If the axes were isotropic, this distribution would be 
uniform. 




Fig. 8. — Median alignment of the major (a-left) and mi- 
nor (b-right) axes at different radii. The alignment is with re- 
spect to the R = 0.06 (black/solid), R = 0.12 (red/dotted), 
-R = 0.25 (green/short-dashed), R = 0.4 (cyan/dot-dashed), R = 
0.6 (blue/dot-dot-dot-dashed), and R = 1.0 rvir (magenta/long- 
dashed) shell. For each pair of radii, only halos where the direc- 
tion of the axis is determined to within 0.2 radians at both radii 
are used. 



of an isotropic error in any particular plane, such as the 
plane containing both of the measured axes, is half of the 
isotropic error, so we divide the isotropic error of 0.2 ra- 
dians by two and add the error of each axis in quadrature 
to find the typical ^spurious ~ 0.14. Therefore, the cosine 
of the angle between the two axes, which is intrinsically 
1.0, is measured to be cos 6'spurious = 0.99. If the axes 
are intrinsically perpendicular, in which case the effect 
of ^spurious on the direction cosine is maximized, the er- 
ror in the direction cosine is 0.14. Most halos have well 
aligned axes (see below), so the error is negligible. 

Figure [3 shows histograms of the alignment between 
the major axis of the R = 0.4 rvir shell and the ma- 
jor axis of the outer (top two panels) and inner (bottom 
three panels) regions of the halo. The alignment is very 
good at all radii. The relative alignment of the major and 
minor axes as a function of radius is shown in Figure |H1 
Each line shows the median alignment with respect to a 
different fiducial radius, recognizable as the radius where 
the median is exactly unity. The axes are extremely well 
aligned within 0.6 rvir- Near the virial radius, some ha- 



8 



Bailin and Steinmetz 




Fig. 9. — Median alignment of the major axis, as in FigurelHk. 
but only for halos with masses of lO" - lO^^ /j-l (if.ftj^ lo"- 
10" h-'^ Mq (middle), and lO" - 3 X lO" h''^ Mq (right). The 
alignment is with respect to the R = 0.06 (black/solid), R = 0.12 
(red/dotted), R = 0.25 (green/short-dashed), R = 0.4 (cyan/dot- 
dashed), _R = 0.6 (blue/dot-dot-dot-dashed), and R = 1.0 rvir 
(magenta/long-dashed) shell. For each pair of radii, only halos 
where the direction of the axis is determined to within 0.2 radians 
at both radii are used. 




Fig. 10. — Histogram of spin parameters A for the halos in our 
sample. The thick line shows the distribution when all of the mass 
in the halo is used to calculate A, while the thin line shows the 
distribution when only the central 0.12 rvir of the halo is used. 
The dashed line is the lognormal fit to the distribution, where 
Ao = 0.035 and cr = 0.58. 



los show deviations, but there is stiU usually very good 
alignment, especially for the minor axis. The alignment 
is better than that seen by JS02. This confirms their sug- 
gestion that many of the halos in which they measured 
poor alignment were nearly prolate or oblate. Such halos 
have large errors in their direction determination, and so 
are not included in our sample. 

Figure 1^1 examines the internal alignment of the major 
axis as a function of halo mass. The inner 0.4 rvir of 
the halos are equally well aligned for halos of all masses. 
However, the outer half of the halo is better aligned with 
the rest of the halo in high mass halos than in low mass 
halos. JS02; saw a similar effect and suggested that it 
was because the low mass halos are intrinsically rounder 
and therefore have larger errors. We rule out this expla- 
nation, as halos with large errors are not included in our 
sample for any mass. Therefore, the stronger alignment 
within high mass halos appears to be a real effect. 

4.2. Angular momentum 

Before we examine the orientation of the angular mo- 
mentum, it is worth commenting on its magnitude in 
our simulated halos. The angular momentum is usually 
quantified by the spin parameter A, where 



l)Peebleslll969() . We use the computationally simpler A' 
as an estimate for A, where 



A' = 



J 



(7) 



V2MVR 

(jBullock et al.ll2001aj) . In Figure EHl we plot the distri- 
bution of spin parameters for our halos. The thick line 
denotes the distribution we find when we use all of the 
mass inside rvir to calculate A. The spin parameters are 
well described by a lognormal distribution, shown as the 
dashed line: 



P(A) = 



1 



AVZttct 



exp 



ln'(A/Ao) 



(8) 
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Fig. 11. — Histograms of the direction cosine between the angular 
momentum of the halo at R = 0.4 rvir a^nd the angular momentum 
at R = 1.0, 0.6, 0.25, 0.12, and 0.06 rvir (top to bottom). Each 
histogram contains halos whose angular momentum direction is 
determined to within 0.4 radians at both radii of the comparison. If 
the orientations were isotropic, this distribution would be uniform. 



with Ao = 0.035 and a = 0.58, in good agreement wit h 
the results of other simulations ijBullock et al.ll2'0 01a?l. 
The distribution of spin parameters remains essentially 
unchanged if we only include mass within the central 
0.12 rvir of the halo, as seen by the thin histogram in 
Figure El indicating that the angular momentum is dis- 
tributed evenly throughout the halo. In the very inner- 
most 0.06 rvir, the spin parameters rise, with the median 
value of A increasing from Aq = 0.035 to 0.047. How- 
ever, we do not believe that this is a physical effect, as 
the particles in this region are affected by the numerical 
force softening, which results in more tangential orbits 
and correspondingly larger spin parameters. 

The orientation of the angular momentum cannot be 
determined as precisely as the orientation of the principal 
axes. In order to have a reasonably large sample of halos. 
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TABLE 2 

Number of 
halos with 

ANGULAR 
MOMENTUM 
DIRECTIONS 
DETERMINED TO 
WITHIN 
0.4 RADIANS 



Radius 


Halos 


0.06 


702 


0.12 


1686 


0.25 


2820 


0.4 


3060 


0.6 


3094 


1.0 


3229 



1.00 7 




0.60 L ^ ^ ^ ^ \ ^ ^ ^ ^ ^ ^ ^ ^ L 

D. I 1 .0 

Fig. 12. — Median alignment of the angular momentum vector 
at different radii. The alignment is with respect to the R = 0.06 
(black/soUd), R = 0.12 (red/dotted), R = 0.25 (green/short- 
dashed), R = 0.4 (cyan/dot-dashed), R = 0.6 (blue/dot-dot-dot- 
dashed), and R = 1.0 rvir (magenta/long-dashed) shell. For each 
pair of radii, only halos where the direction of the angular momen- 
tum vector is determined to within 0.4 radians at both radii are 
used. 



we use angular momentum vectors whose bootstrap error 
is less than 0.4 radians. This is twice as large as the limit 
adopted for the principal axis directions. The number of 
halos that satisfy this criterion at each radius is given in 
Table 121 Following the logic of 14.11 the error in the di- 
rection cosine of two vectors with errors of 0.4 radians is 
0.04 if they are perfectly aligned and 0.28 if they are per- 
pendicular. The median errors of the samples are half of 
these worst-case scenarios. The alignment is shown to be 
good (see below), so the effect of the errors is negligible. 

Figure^]shows histograms of the relative alignment of 
the angular momentum of the R = 0.4 r^iv shell with the 
outer regions of the halo (upper two panels) and the inner 
regions of the halo (lower three panels). The alignment 
is very good at most radii. 

The relative alignment of the angular momentum as 
a function of radius is shown in Figure 1121 Each line 
shows the median alignment of the angular momentum 
with respect to a different fiducial radius, recognizable 
as the radius where the median is exactly unity. The 
alignment gets progressively worse as the radii get fur- 
ther separated; the median cosine between the angular 
momenta in the innermost and outermost regions is 0.64. 




0.1 1.0 0.1 i.O 0.1 1.0 

r/r.,. r/r.„ r/r.„ 

Fig. 13. — Median angular momentum alignment as in Figure ll2l 
but only for halos with masses of lO" - 10^^ h'^ Mq (left), 10~ 
10^3 h-l Mq (middle), and lO" - 3 x lO^-* h-^ Mq (right). The 
alignment is with respect to the R = 0.06 (black/solid), R = 0.12 
(red/dotted), R = 0.25 (green/short-dashed), R = 0.4 (cyan/dot- 
dashed), R = 0.6 (blue/dot-dot-dot-dashed), and R = 1.0 Tvir 
(magenta/long-dashed) shell. For each pair of radii, only halos 
where the direction of the angular momentum vector is determined 
to within 0.4 radians at both radii are used. 




Fig. 14. — Histograms of the alignment between the angular 
momentum vector and the major axis of the R = 1.0, 0.6, 0.4, 
0.25, 0.12, and 0.06 Tvir (top to bottom) shell of each halo where 
the error in the major axis direction is less than 0.2 radians and the 
error in the angular momentum direction is less than 0.4 radians 
at that radius. Due to the symmetry of the axes, the direction 
cosine is always positive. If the orientations were random, this 
distribution would be uniform. 



However, the angular momentum vector at intermediate 
radius, such as at 0.4 Tvir, is generally representative of 
its direction at all radii. 

Figure^] examines the internal alignment of the angu- 
lar momentum as a function of halo mass. The patterns 
seen in Figure ^1 generally hold for all masses, although 
the alignment between the very innermost and outermost 
regions is slightly worse for the highest mass halos. 

4.3. Alignment between the angular momentum and the 

halo shape 

We investigate here the alignment between the angular 
momentum of a halo at a given radius and the principal 
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r/r.» 

Fig. 17. — Median alignment between the angular momentum 
vector and the major (black/solid), intermediate (red/dotted), and 
minor (blue/dashed) axis of each halo as a function of radius within 
the halo. 
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Fig. 15. — As in Figure ICT but for the intermediate axis. 
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Fig. 16. — As in Figures \lAl and [THl but for the minor axis. 



axes of the mass distribution at that radius. If we com- 
pare an angular momentum vector whose error is 0.4 ra- 
dians with a principal axis whose error is 0.2 radians, the 
error in the direction cosine is 0.02 if they are perfectly 
aligned and 0.22 if they are perpendicular. The median 
errors are half of those values. Therefore, in the cases 
where the alignment is good, the effect of the errors is 
negligible. In the cases where the vectors are perpendic- 
ular or isotropic, we must take the errors into account 
when drawing conclusions. 

Figures El and El show histograms of the cosine 
between the direction of the angular momentum vector 
and the major, intermediate, and minor axis respectively. 
Only those halos with both angular momentum direction 
errors of less than 0.4 radians and axis direction errors of 
less than 0.2 radians are used. The angular momentum 
vector tends to be perpendicular to the major and inter- 





FlG. 18. — As in Figure ITtI but only for halos with masses of 
10" - 10^2 ^-1 (left), 10^2 _ iq13 f^-i (middle), and 
10^^ - 3 X 10" h-^ Mq (right). 



mediate axes, and parallel with the minor axis. Because 
of the different effects of the error on parallel and perpen- 
dicular vectors, the tendency of the angular momentum 
to be perpendicular to the major axis is as significant as 
the trend for it to be parallel to the minor axis, despite 
the different appearance of the histograms. The angular 
momentum tends to lie perpendicular to the intermediate 
axis, but this trend is weaker at larger radii. 

These results a r e co nsi stent wi th tho se of 
Barnes fc Efstathioul lll98^ . iDubinskil (|1992l^ . and 
W92L Of these, only IW92I quantify any change with 
radius; they found slightly better alignment at larger 
radii, in contrast to the results presented here. However, 
both of the radii at which they performed the compar- 
ison were well within the virial radius, well inside the 
radii where we see the alignment drop. 

These relationships are summarized in Figure 1171 
which shows the median alignment between the angular 
momentum vector and each of the principal axes. The 
angular momentum tends to lie parallel with the minor 
axis and perpendicular to both the major and interme- 
diate axes. These trends are strongest in the central 
0.25 Tvir of the halo, deteriorating slightly in the outer 
regions. The median cosine of 0.9 corresponds to a mis- 
alignment between the angular momentum and minor 
axis of 25°. Figure [T51 shows how these trends depend 
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Fig. 19. — Projection in xy of the major (top-left), intermediate 
(top-right), and minor (bottom-left) axes, and unit angular mo- 
mentum vectors (bottom-right) of high mass (M > 10^^ h~^ 
halos in a slab of depth Az = 12.5 h~^ Mpc (one quarter of the 
simulation volume). The three-dimensional length of each line or 
vector is 2 h~^ Mpc, including the unseen z component. Dots show 
the positions of all halos, regardless of mass. 



on the halo mass. For high mass halos, the alignment is 
slightly worse and has less of a dependence on the loca- 
tion within the halo. 

5. EXTERNAL ALIGNMENT 
5.1. Introduction 

In this section, we compare the orientation of the prin- 
cipal axes and angular momenta of individual halos with 
the location of mass around them and the orientation of 
those properties in surrounding halos. In §0|we demon- 
strated that the properties of halos in the R = 0.4 Tvir 
shell are characteristic of their properties over a large 
range of radii. Therefore, we use the R — 0.4 rvir mea- 
surements when comparing to other halos. 

For each halo, the volume is split into 7 radial 
bins. The nearest bin spans separations from to 
390.625 kpc. The outer radii double for each sub- 
sequent bin, while the inner radius is equal to the outer 
radius of the interior bin. The largest bin has an outer 
radius of 25 h^^ Mpc, extending to the edge of the pe- 
riodic box. The nominal radius of each bin is the outer 
radius divided by \/2; this corresponds to the geometric 
mean between its inner and outer radius for all but the 
innermost bin. The inner bin has no formal inner radius, 
but in practice is limited by twice the radial extent of 
the typical halo, or 250 kpc. Halos whose centers of 
mass are closer to each other than this are merging, and 
are detected as a single object by the group find e r. 

We follow the nomenclature of iSolinter et alJ l)1997f) . 
and use "alignment" to refer to the tendency of a vector 
(such as a principal axis or an angular momentum vector) 
to point toward or away from other halos, and "correla- 
tion" to refer to the tendency of vectors in different halos 
to point in the same direction. 

A visual impression of how the axes and angular mo- 
menta align can be seen in Figure 1191 which shows the 



axes and angular momentum vectors of the high mass ha- 
los (M > 10^3 h^^ Mq), with the location of the other 
halos shown as dots. Because of the filamentary nature of 
the large scale structure (e.g. fColbcrg ct al. 1999), pos- 
itive alignment indicates that a quantity tends to point 
along filaments. 

5.2. Axis alignments 

We compare here the alignment of the principal axes 
of the halos with the location of surrounding structure. 
Methods of measuring these alignments vary in the lit- 
erature, as does the nomenclature for a given metric. 
We adopt an internally-consistent nomenclature ^^y and 
^\xy\i which are defined to be the mean value of the di- 
rection cosine between direction x and direction y, and 
the mean of the absolute value of the direction cosine re- 
spectively. In general, ^xy is positive if x and y point in 
the same direction and negative if they point in opposite 
directions (this is not applicable if x or y is symmetric, 
such as if it is a principal axis), while ^\xy\ is greater than 
1/2 if a; and y lie parallel to each other and less than 1/2 
if they lie perpendicular to each other. We note names 
other authors use for the same quantities when applica- 
ble. For example, to measure the alignment of the major 
axis, whose direction is defined by the unit vector a, with 
the large scale structure, we calculate the alignment ^\ar\ ■ 

^\ar\{r) = (|a-f|) = ^^|ai -fyl, (9) 

id 

where the sum over i is over all halos in the primary 
sample, the sum over j is over all halos in the secondary 
sample, is a unit vector in the direction of the dis- 
placement from halo i to halo j , and N is the number of 
terms in the double sum. The primary sample consists of 
all halos whose major axis is determined to within 0.2 ra- 
dians, while the secondary sample consists of all halos. 
We define the alignments and ^|cr| for the intermedi- 
ate axis Id and the minor axis c similarly. Note that the 
primary samples used to define ^\ar\j C\br\7 and ^i^^i are 
not identical, as the set of halos with good major axis 
determinations is deficient in very oblate halos, the set 
of halos with good minor axis determinations is deficient 
in very prolate halos, and the set of halos with good in- 
termediate axis determinations is deficient in both very 
prolate and very oblate halos. 

Figure|2nishows histograms of the distribution of direc- 
tion cosines |a • f | for halos at a variety of separations. If 
the axes were randomly oriented, the distributions would 
be uniform. The distribution is mostly isotropic, but 
there is an excess of halos with |a • f | > 1/2, ie. halos 
whose major axes lie parallel to the filaments. We quan- 
tify this by calculating the mean, ^i^j^i (note that this is 
simila r but not identical to the quantity w{r) used by 
IKE04'). 

The upper-left panel of Figure [Til shows ^i^^i for the 
halos. We find that the alignment is well fit by a power 
law over a wide range of separations. The fit is shown as 
the solid line in Figure [21 and is given by 

^\ar\{r) = l+m,r", (10) 

where r is the separation in units of h^^ Mpc, the align- 
ment at 1 h^^ Mpc is mi = 0.015, and the slope is 
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Fig. 20. — Histograms of the direction cosine |a ■ r| between the 
major axes of halos in the primary sample and the displacement 
from each primary halo to all surrounding halos, binned by radial 
separation between the centers of the halos. The direction cosine is 
always positive due to the symmetry of the axes. If the axes were 
randomly oriented, the distributions would be uniform. The radial 
bins consist of halos separated by (top to bottom): 12.5 - 25, 6.25 
- 12.5, 3.125 - 6.25, 1.5625 - 3.125, 0.78125 - 1.5625, 0.390625 - 
0.78125, and - 0.390625 h'^ Mpc respectively. 



0-55 
0-50 




0,55 
50 ■ 



■eporation {ft-' Mpc^ 
i A„ halos 



Separation (/i ' Mpc} 

Fig. 21. — Mean alignment of the principal axes of halos in 
the primary sample with the displacement from the primary halos 
to all surrounding halos, as a function of distance between the 
halos. The alignment of the major axis (5|ar|)i intermediate axis 
(C|6r|)i minor axis (5|cr|) is plotted in the top- left, top-right, 
and bottom-left panels respectively. Only halos where the relevant 
axis is determined to within 0.2 radians are used for the primary 
sample. The different symbols are for primary samples consisting of 
all such halos (crosses) or of only those with masses greater than 
lO'^'^ Mq (diamonds). The diamonds are shifted slightly to 

the right for clarity. Error bars represent the Icr Poisson sampling 
error in the mean. The solid lines are t he p o wer law fi ts fo r the 
full samples, and are given by equations llOi . illl . and 1121 . The 
dotted line is the value expected for random orientations. 



a = —0.64. lF02l find eve n stron ger alignment for their 
cluster-mass halos, as do IKE04I and HBB05| (except at 
the very smallest separations, which is at the spatial limit 
of their simulations and the point where halos begin to 
overlap). This may be because all of these authors ana- 
lyze simulations performed in much larger boxes, which 
contain power on longer wavelengths than exist in our 



smaller simulation volume (which also allows them to 
measure alignments at separations far exceeding our en- 
tire box length). Alternatively, the stronger alignment 
they detect may be a consequence of the higher mass of 
their halos. To see whether the halo mass is important, 
we recalculate C|ar|i restricting the primary sample to 
those halos with masses greater than 10^^ M©. This 
is the tendency for the major axes of group or cluster 
mass halos to point along the filaments. The results are 
shown as the diamonds in Figure !^ The group and clus- 
ter mass halos are much more strongly aligned than the 
full sample (which is dominated by galaxy mass halos). 
The cluster alignment ^i^^i is constant out to 3 Mpc, 
after which it drops until it agrees with the results of the 
full sample by 6 Mpc. HBB05 also find a strong 
mass dependence; however, the alignments of their low- 
est mass sample are still larger than the alignments of 
our highest mass sample, indicating that the box size 
also plays a role. 

The intermediate axes are also aligned with the fil- 
aments, as seen in the upper-right panel of Figure 1211 
While the innermost bins show no alignment, there is 
clear evidence of alignment at separations of 1 Mpc 
and beyond. The solid line is a power law fit to the outer 
5 points, i.e. for r > 781 h^^ kpc, and is given by 



£.\br\ = ^ +mir", 



(11) 



where nii — 0.008 and a = —0.7. The high mass halos 
again show constant alignment out to several Mpc. 

As both the major and intermediate axes tend to point 
along filaments, the minor axis must tend to lie perpen- 
dicular to the filaments. This trend can be clearly seen 
in the bottom-left panels of Figure ^] as sequences of 
parallel lines crossing the filaments. The lower panel of 
Figure 121 shows £,\cr\- The solid line is a power law fit 
for the full sample excluding the innermost bin, i.e. for 
r > 391 kpc, and is given by 



(12) 



where the alignment at 1 Mpc is irii ~ —0.027 and 
a = —0.7. The minor axes of group and cluster mass 
halos show even stronger alignment than do the galaxy 
mass halos, as seen by the diamonds in Figure I2T1 

5.3. Axis correlations 

We compare here the tendency for the principal axes of 
neighbouring halos to point in the same direction. The 
procedure we use is completely analogous to that used 
to calculate the alignments in § 15.21 The correlation be- 
tween the major axes is defined as 



^\aa\ir) = (|a-a|) 



-E 

1,1 



a, ■ ai 



(13) 



where the sum is over all unique pairs (i,j) because the 
measurement is symmetric with respect to each pair of 
halos. Only halos whose major axes are determined to 
within 0.2 radians are used for both the primary and 
secondary samples. 

The top-left panel of Figure |^ shows the mean corre- 
lation ^\aa\ of the major axes as a function of the halo 
separation (this is identical to the quantity defined as 



Alignment of ACDM Halos 



13 



ieporation (ft-' Mpc} 
J This work 

^ Heovens et al. (2000) 

Power low fit 

■■■■ Random 



Seporotion {h ' Mpc} 

Fig. 22. — Mean correlation of the principal axes of halos for 
which the direction of the major axes is determined to within 0.2 ra- 
dians for both halos. The correlations of the major axes (5|aa|)i 
intermediate axes and minor axes (5|cc|) are plotted in the 

upper-left, upper-right, and bottom-left panels respectively. The 
solid line is the po wer law fit for the major axis correlation, and is 
given by equation I14i . The dotted line is the expected value for 
random orientations. Error bars represent the la Pois son samp ling 
error in the mean. Diamonds represent the results of IHRHOCl 



in lKE04j) . Although almost all bins are individually 
consistent with isotropy, they all lie above 0.5, and taken 
together are evidence that the directions of the major 
axes are correlated. We fit a power law and find 

1 

2 



^\aa\{r) 



mir 



(14) 



where r is the separation in units of h"^ Mpc, the cor- 
relation at 1 Mpc is nii = 0.004, and the slope is 
a = —1.2, although the large errors introduce consider- 
able uncertainties in these val ues. The correlations found 
by |F02, KE0| and IHBB05I are considerably stronger 
than those found here. Comparison with the results of 
ICroft & Mctzlcr ( 2000) and HRHOO is more difficult be- 
cause these authors measure projected ellipticity correla- 
tions rather than axis correlations. The ellipticity corre- 
lations are diluted with respect to the axis correlations, 
and therefore we should consider their results as lower 
limits. The results of HRHOO are plotted as diamonds in 
Figure El where we note that 

1 

^\aa\ ^ ^ + — {{eiei) + (6262)) . (15) 

The agreement is excellent. However, because their 
results should be co nsidered as lower limits, and 
ICroft fc Met^ (1200(1 found even larger ellipticity cor- 
relations in the same simulation, it appears that the 
shapes of the halos studied by these authors also cor- 
relate more strongly than the halos in our simulation. 
This is likely because our sample is dominated by galaxy 
mass halos, while all of the previous studies have ana- 
lyzed only cluster mass halos. We have demonstrated in 
S I5.2l that the higher mass halos are more strongly aligned 
with the large scale structure; it is therefore likely that 
they also show stronger correlations with each other. To 
test this, we have recalculated ^\aa\ using only halos with 
masses greater than 10^'^ Mq to see if the behaviour 
of high mass halos differs from those of lower mass, but 
due to the small number of very high mass halos in our 
sample, the errors are too large to draw any conclusions. 
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Fig. 23. — Mean vector alignment between the angular mo- 
mentum vector of a halo and the location of surrounding halos, 
as a function of radial separation. Error bars represent the la 
Poisson sampling error in the mean. The dotted line is the ex- 
pected value for random orientations. Diamonds are shifted to the 
right for clarity. (Top-left): The alignment of all halos (crosses), 
and of only those halos with masses greater than lO"'^ Mq 
(diamonds). (Top-right): The alignment of halos that have 3 or 
fewer neighbouring halos within 2 Mpc (crosses) compared to 
those halos with 4 or more neighbouring halos within 2 Mpc 
(diamonds). The solid line is the power law fit for the low den- 
sity sample at r > 781 h^^ kpc, and is given by equation I17i 
with mi = 0.01 and a = —0.6. (Bottom-left): The alignment of 
halos with concentration parameter Cvir > 16 (crosses) compared 
to those with NFW concentration parameter Cvir < 7 (diamonds). 
The solid line is the power law fit for the high c onc entration sample 
at r > 391 kpc, and is given by equation I17i with mi = 0.03 
and a = —1.3. 



The mean correlations of the intermediate and minor 
axes of the halos, ^i^^i and ^|cc|> ^re plotted in the upper- 
right and lower- left panels of Figure [221 The error bars 
are too large to robustly detect any correlation, though 
the lower panel shows a suggestive preponderance of bins 
with ^|cc| > 0.5. A larger sample of halos simulated at 
equally high spatial resolution is required to confirm if 
this is real. 

5.4. Angular momentum 

We investigate the tendency for the angular momen- 
tum to point toward or away from nearby halos, for it to 
lie parallel or perpendicular to the filamentary structure, 
and for the angular momentum vectors of neighbouring 
halos to point in the same direction. To measure the ten- 
dency of halo angular momenta to point toward or away 
from other halos, we calculate the alignment 



(16) 



where Li is a unit vector in the direction of the angu- 
lar momentum vector for halo i. The primary sample 
consists of all halos where the direction of the angular 
momentum is determined to within 0.4 radians, while 
the secondary sample consists of all halos. We show the 
results in the top-left panel of Figure [23 There may be 
a weak tendency for the angular momentum to point 
toward local density enhancements on scales of 0.5 - 
2> Mpc, but the size of the error bars makes such 
a result tentative. When we restrict the sample to those 
halos with masses greater than 10^^ M©, we find no 
statistical deviation from random. 
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iKashikawa &: Okamural l)1992|) and iNavarro et al.l 
11200411 have found that galaxies within the local super- 
cluster have their spins pointing within the supercluster 
plane, while those at least 2 Mpc from the plane have 
spins that point toward or away from the plane. To sec if 
the behaviour of simulated halos in low and high density 
regions differs, we have split the sample into those halos 
that have 3 or fewer neighbouring halos within 2 Mpc 
(the low density sample), and those that have 4 or more 
neighbours within 2 Mpc (the high density sample). 
There are 2155 and 1714 halos in the low and high density 
samples respectively. In the top-right panel of Figure ESI 
we have plotted £,Lr for these two samples. By construc- 
tion, there are very few pairs (and therefore large error 
bars) at small separations in the low density sample, but 
beyond 1 Mpc there is a clear detection of positive 
alignment in this sample. The alignment is well fit by a 
power law of the form 

a,- = mir", (17) 

where r is the separation in units of Mpc, the cor- 
relation at 1 Mpc is mi = 0.01, and the slope is 
a = —0.6. The high density sample shows no coherent 
tendency for the angular momentum vectors to point to- 
ward or away from density enhancements. The low and 
high density samples do not show significant deviations 
from the full sample for any of the other statistics stud- 
ied. 

In order to see how the internal structure of the 
halo affects its alignment properties, we have split the 
sample int o those halos with particularly high NFW 
ijNavarro e t al. 1996j) concentration parameters (cvir > 
16) and those with particularly low concentration pa- 
rameters (cvir < 7). As demonstrated by the bottom- 
left panel of Figure there is a clear signal of pos- 
itive alignment in the high-concentration sample, and 
none in the low-concentration sample. The alignment 
of the high-concentration sample is well fit by a power 
law of the form of equation H17|l . with mi = 0.03 and 
a — —1.3. Because halo concentration anti-correlates 
with mass (Bullock et al. 2001b), our high-mass sample 
is deficient in high-concentration halos. Therefore, the 
tendency for low mass halos to show positive ^Lr while 
high mass halos show no such a tendency is probably 
a direct consequence of the correlation between and 

Cvir ■ 

The tendency for the angular momentum vectors to lie 
parallel versus perpendicular to the filaments is measured 

by 

f|i,|(r)^(|L.f|)^l^|L,.f,,|, (18) 

defined as L|| in lHatton fc NininI 1)200 Ij) . or 

coir) ^ (|L . f p) - i ^ 1 ^ |L, . f,, p - 1, (19) 

as used by ILee fc Peiil 1)200 1)) . We plot these quanti- 
ties in Figure |^ The solid line in the b ottom panel 
shows the prediction of ILee fc Penl l)2001)) from linear 
tidal torque theory. On scales less than 3 h^^ Mpc, 
the angular momentum tends to lie parallel to the fil- 
am ents. The valu e s are consistent with t hose found 
bv iHatton fc NininI l)2001)) . However, both iLee fc Penl 
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Fig. 24. — Mean alignment I (top) and uj{r) (bottom) hetween 
the angular momentum vector of primary halos and the location 
of all surrounding halos. The different symbols are for primary 
samples consisting of all such halos (crosses), or of only those with 
masses greater than 10^^ Mq (diamonds). The diamonds are 
shifted slightly to the right for clarity. Error bars represent the la 
Poisson sampling error in the mean. The dotted line is the expected 
value for random orientati ons. The solid lin e is the prediction from 
linear tidal torque theory ILee fc Penll200ll) . 
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Fig. 25. — Mean correlation ^i^j^ of the angular momentum vec- 
tors of halos as a function of their separation. Error bars represent 
the Icr Poisson sampling error in the mean. The dotted line is the 
expected value for random orientations. 

1)2001)) and IF02I find that the angular momenta of halos 
tend to lie perpendicular to the filaments. This discrep- 
ancy m ay lie in the different mass ranges probed. The 
halos of lHatton fclji nin (2001') cover a very similar mass 
range to those in this work, while the sample of FQ2 con- 
sists entirely of cluster mass halos (their smallest halo 
has a mass of 1.4 x 10"'^'* M©, nearly the mass of our 
largest halo). We have recalculated ^\Lr\ and uj{r) using 
only halos with masses greater than 10^"^ Mq in the 
primary sample and plotted them as diamonds in Fig- 
ure |31 The behaviour of the high mass halos is radically 
different; the angular momenta of groups and clusters 
tend to point perpendicular to the filaments. 

The correlation between the directions of the halo an- 
gular momenta is defined as 

aL^(L.L)^l^L,.L„ (20) 

where the sum is over all unique pairs (i,j). This 
is equivalent to th e quantity defined is r](r) in 
IPorciani et all ll2002l . As seen in Figure we detect 
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Fig. 26. — Mean correlation (top) and r;(r) (bottom) be- 

tween the angular momentum vectors of halos as a function of 
their separation. Error bars represent the la Poisson sampling 
error in the mean. Th e solid line shows the analytic result of 
ICrittenden et aTl 120011) . where we set the correlation parameter 
between the shear and inertia tensors a = 0.55, while the solid 
line is the linear tidal torque prediction of lPen et al.l 120001) . The 
dotted line is the expected value for random orientations. 



no deviations from random. We also measure the ten- 
dency for halo angular momentum vectors to be either 
parallel or perpendicular to each other. We define two 
quantities to measure this, 

^l^i|^(|L.L|)^l^|L,.L,|, (21) 

which is t he same as the quantity defined as /i(r) in 
iHatton fc Ninin ( 2001), and 

rKr) ^ (|L • Lp) - i ^ 1 ^ |L, • L,f - i (22) 

which is the same as iPen et all ||200(TI) 's r]{r) and 
iPorciani et alJ ^QQMi's mif)- The quantities S,\LL\{r) 
and 'q{r) are shown in Figure EEl We detect no devia- 
tions from iso tropy. For referen c e, we also plot the linear 
predictions of lCrittenden et all l)2001[l (where we set the 
correlation parameter between the shear and inertia ten- 
sors a = 0.55 so that the linear results al so match the 
TV-body results of HRHOQl and iPen et al] IpjOQI . Due 
to the size of the error bars, our results are consistent 
with the non-detections and weak correlations fo und by 
other authors IhRHOO: Pen et al. 2000; Crittende n et all 
|2001: Hatton & Ninin 2001t iPorciani et al.ll2002atlF0lK 
The number of halos in our sample is too small to detect 
angular momentum correlations in any particular sub- 
sample. 

6. DISCUSSION 

6.1. Internal alignment and galactic warps 

The internal alignment of the angular momentum pro- 
vides several clues to the forces a galaxy undergoes dur- 
ing its evolution. In particular, we note that, as demon- 
strated in Figure El the angular momentum and the mi- 
nor axis of a halo are typically misaligned by ~ 25° , even 
in the central regions of the halo where a galactic disk 
would lie. Therefore, if the angular momentum of the 
dark matter and baryons are aligned, and the disk lies 
perpendicular to its angular momentum, then the disk 



and halo will be typically misaligned by ~ 25°. Such mis- 
alignment could^_be_2ajtli;_res2onsible for warped galac- 
tic di sks l|Dekel fc ShlosmaiJll983t lToomrelll98l iBaihnI 
120041) . However, simulations of live disks in live halos 
demonstrate that misaligned disk-halo sy stems quickly 
reach an unwarped realigned equilibrium l|Binnev et alJ 

However, if the orientation of the disk and/or the ori- 
entation of the halo chan ges, then such an equilib rium 
may never be established. iQuinn fc BinnevI (|1992) sug- 
gested that the infalling angular momentum may be un- 
correlated with the orientation of the central angular mo- 
mentum. Figure 1121 demonstrates that this is true in 
our halos; while the angular momentum at intermedi- 
ate radius is generally representative of its orientation at 
all radii, the angular momentum of the innermost and 
outermost regions are nearly uncorrelated. Therefore, 
as the disk accretes angular momentum, its orientation 
change s. This reorientation of the disk may itself cause 
warp s l)Ostriker fc BinnevlllQSfll : iLope z-Corre doira et alJ 
120021) , or simply prevent the disk and halo from reaching 
an ali gned equilibrium. Similarly, B ailin fc Steinmetd 
(|200l found that over 90% of dark matter halos show 
figure rotation; this reorientation of the halo also acts 
to maintain misalignment between the disk and halo. 
Therefore, it appears that misalignment between the halo 
and disk is common, and perhaps responsible for galactic 
warps. 

6.2. External alignment 

The alignment of the halo axes with the large scale 
structure seen in Figure [Til provides several clues to the 
origin of halo orientations. For instance, the relative 
alignment of the major, intermediate, and minor axes 
is intriguing. The relative strength of the alignment is 
well described by mi, the value of the power law fit at a 
separation of 1 h'^ Mpc. This is 0.015, 0.008, and -0.027 
for the major, intermediate, and minor axes respectively. 
While all previous authors have neglected the minor axis, 
we find that its alignment perpendicular to filaments is 
stronger than the alignment of the major axis along fil- 
aments. Therefore it is the minor axis, not the major 
axis, that is most influenced by the presence of surround- 
ing material. This also explains how both the major and 
intermediate axes can be positively aligned with the fila- 
ment: if the minor axis of a halo lies perpendicular to the 
filament, then both the major and intermediate axes are 
constrained to lie within a plane that contains the fila- 
ment. Therefore, they are both more likely to point along 
the filament than a randomly-oriented three-dimensional 
axis, and therefore both show positive alignment (for ex- 
ample, if all of the minor axes were perfectly perpendic- 
ular to the filament and there were no difference between 
the major and intermediate axes, then ^|cr-| would van- 
ish and both S_\ar\ and ^\br\ would equal 2/7r w 0.64). 
The geometry of the environment is not strictly linear, 
especially at larger scales. However, despite these com- 
plications, the relationships we have found provide use- 
ful quantitative predictions that can be tested when full 
three-dimensional observations of halo shapes are avail- 
able. 

Previous studies have found that the orientation of the 
major axis of a cluster is strongly affected by the di- 
rection of the most recently accreted subhalo, as seen 
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in simulation s ijvan Haarlem fc van de Wevgaerti Il993t 
lTormenllT997|) and observations l|Ebeling et al.ll2004^ . If 

these are accreted from filaments, then there should be 
a strong alignment of cluster major axes with the fil- 
aments, as we have found. The minor axes are even 
more strongly aligned; this either suggests that recently 
accreted subhalos have a stronger effect on the minor 
axis than on the major axis, or more likely that figure 
rotati on scrambles this alignment. Bailin & Stcinmet^ 
(|2p04') found that most halos (78%) show slow figure ro- 
tation about the minor axis, with a median pattern speed 
of w 0.15 h km s~^ kpc~^, suggesting that the major axis 
of a typical halo can change by 90° in a Hubble time. 
However, only 13% of halos show figure rotation about 
their major axis. Therefore, even if a halo is formed with 
both its major and minor axes perfectly aligned with the 
filament, the major axis most often rotates away from 
that orientation by several degrees while the minor axis 
most often remains in its original orientation. 

The galaxy mass halos show weaker alignment between 
their axes and the large scale structure than do the 
cluster mass halos. It may be that the axes of lower 
mass halos are less affected by recent major mergers, 
or that the direction of the accretion i s mor e isotropic 
l)Vitvitska et alJl2002() . lAubert et all 1)20041) find that 
the accretion onto halos at masses down to 5 x 10^^ il/© 
is quite anisotropic; however, 97% of our "galaxy" mass 
halos lie below this limit. On the other h and, figure rota- 
tion m ay again be at work. Although Bail in fc Steinme't^ 
l|200l found no relationship between the figure rotation 
pattern speed and the halo mass, lower mass halos are on 
average dynamically older, and therefore figure rotation 
may have had a longer time to modify the original ori- 
entation. In any case, although we cannot directly test 
whether the correlation of the axes depends on mass, 
the strong dependence of axis alignment on mass sug- 
gests that the intrinsic alignments of galaxy mass halos 
are much weaker than the intrinsic alignments of cluster 
mass halos. Therefore, when interpreting weak lensing 
results, it is important to note that the predictions for the 
intrinsic alignments of galaxies ma de by analyzing c lus- 
ter m ass halos fe.g. lCroft fc Metzle3l2000HHeavens et al.l 
l200nt) are likely overestimates. 

Turning to the angular momentum, we find an intrigu- 
ing correlation between the concentration of a halo Cvir 
and the tendency for its angular momentum to point 
toward other halos, ^Lr (see Figure I23II . While most 
halos show no tendency for the angular momentum to 
point toward rather than away from nearby halos, ha- 
los with particularly high concentration parameters show 
positive ^Lr- The concentration parameter is a measure 
of the fa ctor by which the halo has collapsed since its for- 
mation ijNavarro et al. 1997) ; therefore non- linear effects 
and higher order terms in the expansion of the poten- 
tial are expected to play a more important role in these 
halos. Indeed, in simple linear theory where the poten- 
tial is approximated by a second-order Taylor expansion 
ifWhite .198^ . each axis is symmetric and so it is im- 
possible to produce a preferred direction. However, a 
second-order expansion of the potential $ is only a solu- 
tion to the Poission equation in a constant density back- 
ground. In order for the angular momentum to point 
toward nearby structure, there must be a background 
density gradient. In the simplest case, the potential must 



be expanded to third order to account for a linear den- 
sity gradient. Each additional derivative of the potential 
couples to additio nal moment s of the density distribution 
of the protohalo l|Porciani et al. 2002a). Therefore, the 
third derivatives of the potential couple to the skewness 
of the density distribution. Because the third derivatives 
are no longer symmetric (in general Taj ^ Tijj, where 
Tijk = d^^/dqidqjdqk), the symmetry is broken and it 
is possible to produce a preferred direction. 

The differences between the alignment of the angular 
momentum of low and high mass halos are also intrigu- 
ing. Angular momentum is usually thought to either 
arise from the tidal torquing of an asymmetric protohalo 
|Whitc 1983), or by the accretion of substructure o n non- 
radial orbits l)Tormenlll997t iVitvitska et a l. 2002). In a 
sense, these are not distinct scenarios; accreted subhalos 
are protogalactic material that has been tidally torqued. 
However, the detailed orientation of the angular momen- 
tum may depend on the dumpiness of the acc retion and 
other non- linear effects (Por ciani et al.ll2002al) . 

We find that the angular momentum in galaxy mass 
halos points parallel to the large scale structure, while 
the angular momentum in cluster mass halos points per- 
pendicular to the large scale structure (see Figure I24|l . 
If the direction of the angular momentum is dominated 
by mergers of objects falling in along filaments, then the 
angular momentum will tend to point perpendicular to 
the velocity of the infalling objects, i.e. perpendicular 
to the filament. Therefore, the difference between the 
alignment of the angular momentum of the low and high 
mass halos suggests that the accretion onto low mass 
halos is either smoother than the accretion onto high 
mass halos, or that a significant fraction of the accretion 
in low mass halos occ urs perpendicular to the filament 
ijVitvitska et alJl2002|) . The group and cluster mass ha- 
los, on the other hand, tend to occur at the intersections 
of large filaments, and have been built up by the recent 
accretion of halos along the filament which contribute 
angular momentum perpendicular t o the filamen t (note 
that this confirms the suggestion of iMaller et alJ 1)20021 ) 
that accretion along a preferred axis is required to ex- 
plain the distribution of spin parameters seen in simu- 
lated halos). 

This difference between the inferred origin of the angu- 
lar momentum of low and high mass halos meshes nicely 
with that of Pcirani ct al. (2004). These authors found 
that halos which have experienced a major merger ac- 
quired their angular momentum much more rapidly than 
those halos which grew only through steady accretion, 
suggesting that the sources of their angular momenta 
are the merger events. The halos in the merger sample 
also built up their mass more quickly due to those merg- 
ers, and consequently have larger median present-day 
masses. Our group and cluster mass halos thus contain a 
larger fraction of major merger products, and therefore 
a larger number of halos whose angular momenta were 
built up through mergers. It is encouraging that our re- 
sults, based on the orient ation of the angular momentum 
vectors, and the results of lPeirani et al.lp004() . which are 
based on the evolution of the magnitude of the angular 
momentum, build a consistent picture. 

It is interesting to compare our result s for coUision- 
Icss dark matter halos with the results of iNavarro et alJ 
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l)2004(l for sim ulations that include b aryonic physics. The 
main results of lNavarro et al.l l)2004f l are that the angular 
momenta of baryonic galactic disks in gasdynamical sim- 
ulations tend to align with the intermediate axis of the 
local baryon distribution on scales of « 2 Mpc, and 
that observed edge-on disk galaxies in the local super- 
cluster have their spin axes lying within the supergalac- 
tic plane, as expected in such a scenario. We find that on 
the scale of an individual halo, the angular momentum of 
the dark matter aligns with its minor axis. However, we 
also find that it lies parallel to the large scale structure, 
as expected if the spin ve ctor lies within local sheet-like 
structures as suggested bv lNavarro et alJ ({20041 ^. We also 
find that in low density regions, halo spins point toward 
nearby filaments, and are therefore aligned with the in- 
termediate axis of the local density field. Therefore, de- 
spite the common presence of misalignments between the 
angular momenta of t he baryons and the dark matter in 
individual galaxies ( Sh arma fc Steinmetdl20fl3) . each re- 
tains some memory of the initial torques provided by the 
large scale structure. 

Unlike for galaxy mass halos, for which we can con- 
firm that the angular momenta of the dark matter and 
baryons share similar relationships to the large scale 
structure in simulations, we do not have high resolution 
gasdynamical simulations of the more massive groups 
and clusters with which to compare. It may be possi- 
ble to determine the rotation axis of baryons in these 
systems observationally. In theory, the spin vector of a 
relaxed cluster can be deduced from the presence of a 
redshift gradient of the galaxies in the cluster; however, 
confusion due to structure along the line of sight, the 
small magnitude of the rotation compared to the intrin- 
sic velocity dispersion, ambiguity in the orientation of the 
ellipsoidal shape, intrinsic distance gradients, the lack of 
a large sample of substructureless relaxed clusters, and 
the uncertain relationship between the angular momen- 
tum of the cluster galaxies and that of the smooth X-ray 
emitting gas that dominates the baryonic mass make this 
measurement difficult. The rotation of the X-ray gas 
itself, however, may be measured using the kinematic 
Sunyaev-Zeldovich effect in future Cosmic Microwave 
Background (CMB) s urveys (jChhiba, Ma,rmheim..2nn2; 
ICooravfc Che;il2nnl. 

7. SUMMARY 

We have studied the internal shapes and angular mo- 
menta of galaxy and group mass dark matter halos 
formed in a ACDM A^-body simulation, and studied how 
they are correlated with the large scale structure and the 
properties of neighbouring halos. 

Internally, halos are triaxial with b/a and c/a ratios of 
0.75 ±0.15 and 0.6 ±0.1 respectively. The distribution of 
axis ratios has a tail to low values. The two-dimensional 
projected ellipticities cover a broad range of values from 
to 0.5, with a mean of 0.24, consiste nt with the weak 
lensing results of ' Hoekstra et all l)2004j) . The axis ratios 
rise between 0.12 and 0.6 rvir, beyond which they drop. 
Within 0.12 rvir, the measurement is probably compro- 
mised by the force softening in the simulations. Halos 
are most often prolate in the inner regions, but tend to 
a more even mix of prolate and oblate at large radii. 

The internal alignment of the halos within 0.6 rvir is 
very good, particularly for the minor axis, with a slight 



decrease in alignment in the outermost regions. High 
mass halos have particularly well-aligned axes. While 
the orientation of the angular momentum is also rela- 
tively constant, it changes more noticeably as a function 
of radius than do the axes; there is very little correla- 
tion between the very innermost regions and the very 
outermost regions. The distribution of spin parameters 
is essentially independent of the size of the region used 
to calculate it, indicating that the relative importance of 
the angular momentum is the same at all radii. At any 
given radius, the angular momentum vector tends to be 
aligned with the minor axis and be perpendicular to the 
major axis. However, the mean misalignment of 25° 
implies that galactic disks are generally misaligned with 
their halo. Figure rotation of the halo and reorienta- 
tion of the disk due to newly accreted angular momen- 
tum act to maintain this misalignment. The alignment 
between the angular momentum and the principal axes 
gets weaker at larger radii. The properties of the halo 
at 0.4 rvir SLie quite characteristic of their values at most 
other radii. 

The minor axes of halos show a strong tendency to 
lie perpendicular to the filaments. As a consequence, the 
major (and, to a lesser degree, intermediate) axes tend to 
point along filaments. Figure rotation about the minor 
axis may be responsible for the smaller degree of major 
axis alignment. These alignments fall off with distance 
as a power law. In all cases, the alignment for group 
and cluster mass halos is stronger and extends to much 
larger separations than for galaxy mass halos. Therefore, 
previous studies which predict the intrinsic alignments of 
galaxies based on the shapes of cluster mass halos appear 
to systematically overestimate the correlations. The ma- 
jor axes of halos show a weak correlation with those of 
other nearby halos, while there is no robust detection of 
a correlation for the intermediate or minor axes. 

The angular momenta of high-concentration halos tend 
to point toward, rather than away from, nearby ha- 
los. These halos have collapsed further into the non- 
linear regime, and therefore the effects of higher-order 
derivatives of the tidal field, which are necessary to 
create a preferential direction, are more important for 
these halos. The angular momenta of halos in low 
density environments also tend to point toward local 
density enhancements, in agreement with the results of 
^ishikawa & Okamura (1992) for galaxies 2 Mpc or 
more away from the local supergalactic plane. The angu- 
lar momenta of galaxy mass halos show a weak tendency 
to point along filaments on scales up to 3 Mpc, but 
those of group and cluster mass halos show a very strong 
tendency to point perpendicular to the filaments. This 
appears to be due to the different merger histories of the 
two samples; higher mass halos are more likely to have 
experienced a major merger along a filament which dom- 
inated the evolution of their angular momentum, while 
lower mass halos are more likely to have had a smoother 
accretion history. We detect no correlations between the 
angular momentum directions of nearby halos, but due 
to the size of the error bars, this is consistent with pre- 
vious linear and A'-body studies that predict weak cor- 
relations. Comparisons with recent gasdynamical simu- 
lations and observations of edge-on disk galaxies in the 
local supercluster suggest that both the baryons and dark 
matter in galaxies share a memory of the orientation of 
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the large scale structure that provided the initial torque. 
The alignments we predict may be tested with large sam- 
ples of galaxy redshifts within relaxed clusters, or by 
kinematic Sunyaev-Zeldovich studies in future CMB ex- 
periments. 
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